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O ■ 

(^ ■ Abstract. Wc construct an explicit solution of the Cauchy initial value problem for certain diffusion- 

CNJ I type equations with variable coefRcients on the entire real line. The corresponding Green function 

hjQ' (heat kernel) is given in terms of elementary functions and certain integrals involving a characteris- 

^ I tic function, which should be found as an analytic or numerical solution of the second order linear 

.^^ ■ differential equation with time-dependent coefficients. Some special and limiting cases are outlined. 

^^ I Solution of the corresponding non-homogeneous equation is also found. 

S: 

r-| ! 1. Introduction 

^ ■ In this paper we discuss explicit solution of the Cauchy initial value problem for the one- 

'— '. dimensional heat equation on the entire real line 

(^ . where the right hand side is a quadratic form Q {p, x) of the coordinate x and the operator of 
"^ ■ differentiation p = d/dx with time-dependent coefficients; see equation (12 .ip below. The case of a 
t^^ ■ corresponding Schrodinger equation is investigated in [B] . In this approach, several exactly solvable 
r^ ■■ models are classified in terms of elementary solutions of a characterization equation given by (12.131) 
O ■ below. Solution of the corresponding non-homogeneous equation is obtained with the help of the 
j> I Duhamel principle. These exactly solvable cases may be of interest in a general treatment of the 
nonlinear evolution equations; see |3], [1], [5], [26] and references therein. Moreover, these explicit 
solutions can also be useful when testing numerical methods of solving the semilinear heat equations 



X 

. 5t ■ with variable coefficients. 

2. Solution of a Cauchy Initial Value Problem: Summary of Results 
The fundamental solution of the diffusion-type equation of the form 

-Q^ = ^it)Q^-Ht)x\ + c{t)x— + dit)u + f{t)xu-g{t)—, (2.1) 

where a{t) , b (t) , c(t) , d(t) , f (t) , and g (t) are given real- valued functions of time t only, can be 
found by a familiar substitution 

u = Ae^ = A (t) e^(^'^'*) (2.2) 
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with 

A = A{t)= , ^ (2.3) 

and 

S = S{x,y,t) = a (t) x^ + (3 (t) xy + 7 (t) y^ + 5 (t) a; + e (t) y + k (t) , (2.4) 

where a{t) , (3 (t) , 7 (t) , (5 (t) , e (t) , and k (t) are differentiable real-valued functions of time t only. 
Indeed, 

- = a(-j -6x^ + /a:+(cx-,)- (2.5) 



-a^^ - d = -2a (t) a (t) - d (t) . (2.6) 



provided 

2/i 9x^ 

Equating the coefficients of all admissible powers of x"^y"' with < m + n < 2, gives the following 
system of ordinary differential equations 

dcy. 

— + 6(t) -2c(t)a-4a(t)a^ = 0, (2.7) 

dS 

^-(c(t)+4a(t)a(t))/3 = 0, (2.8) 

^ - a (t) /^^ (t) = 0, (2.9) 

- - (c (t) + 4a (t) a (t)) <5 = / (t) - 2« (t) <? (t) , (2.10) 

^ + (^(t)-2a(t)5(t))/?(t) = 0, (2.11) 

^ + g(t)6{t)-a{t)6'{t) = 0, (2.12) 

at 



where the first equation is the familiar Riccati nonlinear differential equation; see, for example, |12j . 
[T8] . |22j . [23] . [27] and references therein. 

4aa' + 4a6 - 2c (4aa) - (4aa)^ = 0, 4aa = -2d - — 



from 02. 7p and (12.61) and the substitution 



4aa' = -2d'-^^+(^^]\^(2d+i^ 



results in the second order linear equation 

fi" - r (t) ft' - Aa (t) ft = (2.13) 

with 

r(t) = — + 2c-4:d, a(t) =ab + cd-d^ + -l—-—] . (2.14) 

a 2 \a d J 

As we shall see later, equation (12.131) must be solved subject to the initial data 

//(0) = 0, /i' (0) = 2a (0) ^ (2.15) 
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in order to satisfy the initial condition for the corresponding Green function; see the asymptotic 
formula (12.241) below for a motivation. Then, the Riccati equation (12. 7p can be solved by the back 
substitution (12.61) . 

We shall refer to equation (I2.13P as the characteristic equation and its solution /i (t) , subject to 
(I2.15p . as the characteristic function. As the special case (I2.13P contains the generalized equation 
of hypergeometric type, whose solutions are studied in detail in ^20j; see also [1], [19], [25], and [27] . 

Thus, the Green function (fundamental solution or heat kernel) is explicitly given in terms of the 
characteristic function 

U = K (x V t) = ^ ^a{t)x^+l3(t)^y+'y(t)y2+5{t)x+e{t)y+K(t) (2 16) 



Here 



with 



(,, i_fm_diti 

^^~ 4a(t)/i(t) 2a {ty ^ ^ 

/? (t) = ^ exp n (c (r) - 2d (r)) dA , (2.18) 

,(t) ^ --^ exp (2 /\c(.)- 2. (.))..) (2.19) 



1 



6{t) = — -exp( / (c(r)-2rf(r)) rfrj / exp ( - / {c (X) - 2d (X)) dX 






'j^S it) exp (£ (c (r) - 2d (r)) dr^ (2.21) 

'a{r)a (r) 



Kit) 



..2 exp / (c(A)-2rf(A)) dX]{f,{r)6{r)) dr 
10 (Ai'(r)) \Jo ) 

+2 ^* -^ exp [J^ (c (A) - 2d (A)) rfA^ (^/ (r) + ^^ (r) ) rfr, 
a it) fi (t) 



,,.. SHt)-^t'-}^^{,ir)6ir)fdr (2.22) 

+2 r ^ (/. (r) 5 (r)) f / (r) + ^^g {r)) dr 

Jo /^ M V a{T) J 

HO) = ^y e{0) = -6{0), k{0) = 0. (2.23) 
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We have used integration by parts in order to resolve the singularities of the initial data; see section 3 
for more details. Then the corresponding asymptotic formula is 

as t -^ 0^. Notice that the first term on the right hand side is a familiar heat kernel for the diffusion 
equation with constant coefficients (cf. Eq. (15.21) below). 

By the superposition principle, we obtain solution of the Cauchy initial value problem 

— = Qu, u{x,t)\^^Q = uo{x) (2.25) 

on the infinite interval — oo < x < oo with the general quadratic form Q {p, x) in (12. ip as follows 

u{x,t)= / K{x,y,t) uo{y) dy = Hu{x,0). (2.26) 



This yields solution explicitly in terms of an integral operator H acting on the initial data provided 
that the integral converges and one can interchange differentiation and integration. This integral is 
essentially the Laplace transform. 



In a more general setting, solution of the initial value problem at time to 

du 
'dt'' 
on an infinite interval has the form 



Qu, u{x,t)\^^^^ = u{x,to) (2.27) 



u{x,t)= K{x,y,t,to) Uo{y,to) dy = H{t,to)u{x,to) (2.28) 

J —oo 

with the heat kernel given by 

K (x V t tn) = ^a{t,to)x^+(3{t,to)xy+'Y{t,to)y'^+5(t,to)x+e{t,to)y+K{t,to) ^2 29) 

A/27r/i (t, to) 

The function n{t) = jj, (t, to) is a solution of the characteristic equation (12.131) corresponding to the 
initial data 

/i(to,to)=0, /i'(to,to) = 2a(to)^0. (2.30) 

If {/ii,/i2} is a fundamental solution set of equation (I2.13p . then 

^ ^^' ^°^ " wi!il^\^) ^^' ^^°^ ^' ^^^ ~ ^' ^^^ ^' ^^'^^ ^^'^^^ 

and 

/i' (t, to) = J""^^'^ ifi, (to) /i'2 (t) - fi[ (t) /i2 (to)) , (2.32) 

where W (/i^, /X2) is the value of the Wronskian at the point to. 

Equations (I2.17p - (l2.22p are valid again but with the new characteristic function fi (t, to) . The 
lower limits of integration should be replaced by to. Conditions (I2.23p become 

6{to,to) = -e{to,to) = Pj^, ft:(to,to) = (2.33) 

2a (to) 
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and the asymptotic formula fl2.24p should be modified as follows 

^S{x,y,t,to) 

K{x,y,t,to) = . , . (2.34) 

A/27r/i(t,to) 

1 ( {x-yf \ f g (to; , ^ , 



^/4na (to) (t - to) \ 4a (to) (t - to)y "^ \2a{to) 

We leave the details to the reader. 

3. Derivation of The Heat Kernel 

Here we obtain the above formulas (I2.17I) - (I2.22I) for the heat kernel. The first equation is a direct 
consequence of (12.61) and our equation (12. 8p takes the form 

(/x/3)' = (c-2rf)(/i/5), (3.1) 

whose particular solution is (12.181) . 

From (EH) and fl2T8D one gets 

^ W = / 4^^'"^*^ ^^' ^ W = r (^ (^) - 2^ (^)) ^^ (3-2) 

J f^ iV Jo 



and integrating by parts 



7W = - / ^^ - = r+ (^) -■ 3.3 



But the derivative of the auxiliary function 



^ W ^2h{t) 



F{t) = ^ e''^'^ (3.4) 

is 

p' u) = (^' + 2^'«) e^V - Qe^V ^ 4aa/x ^fe ^ 4a/i^ 

in view of the characteristic equation fl2.13p - (l2.14p . Substitution into (13.31) results in (12.190 . 
Equation (I2.10p can be rewritten as 

{^e^^S)' = fie^^{f-2ag), h= f {c~2d) dr (3.6) 

Jo 

and its direct integration gives (I2.20p . 

We introduce another auxiliary function 

G (t) = /i (t) 5 (t) e-'^W (3.7) 

with the derivative given by (13. 6p . Then equation (12. lip becomes 

dt II fi 

and 

e(t) = - f ^e^ dt + 2 f '^e^^ dt. (3.8) 

J ^^ i /^^ 
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Integrating the second term by parts one gets 



FG f (FG)' 



,.^^"* - -] f^'" %)-]''' %^ P-«) 



dt. 



H J fx 
where 



(FG)' = F'G + FG' (3.10) 

- '^-^Me'^ + ^eV + ^e^ + y 



in view of (13. 5p and (13.61) . Then substitution (I3.10p into (13.91) allows to cancel the divergent integrals. 
As a result one can resolve the singularity and simplify expression (13. 8p to its final form (I2.2ip . 

Finally, by fl2:T2|) and ([3TD 

K{t) = - f g6dt+ f^e^^ dt, (3.11) 

where the last integral can be transformed as follows 

aG^ 2h , [ ^^2 J^\ FG^ f (FG^)' , 

-^-e^'' dt = - FG^ di-] = + / ^ '- dt 3.12 

/^^ J V/^y f^ J f^ 

with 

{FG^y = F'G^ + 2FGG' = (FG)' G + FGG' (3.13) 

4aa/i , ^,2 2a/i 2rf/i 



(/iO ^i' '"^ " /i' 

Substitution (I3.12l) - p.l3l) into (13.111) gives our final expression (I2.22p . 

The details of derivation of the asymptotic formula (I2.24p are left to the reader. 

4. Special Initial Data 

In the case u (x, 0) = Uq =constant, our solution (12.260 takes the form 

u{x,t) = / K{x,y,t) uq dy (4.1) 



J —CO 



Uq 



:exp 



(4a7 - (3^) x^ + 2 (27^ - (3e)x + 47/? - e' 



V^Wn \ 47 

provided 7 (t) < 0, with the help of an elementary integral 

The details of taking the limit t -^ 0^ in (14.11) are left to the reader. 



oo 
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When u {x, 0) = 6 {x — Xq) , where 6 (x) is the Dirac deha function, one gets formally 

u{x,t)= / K{x,y,t) 6{y-Xo) dy = K {x,Xo,t) . (4.3) 



Thus, in general, the heat kernel fl2.16p provides an evolution of this initial data, concentrated 
originally at a point Xq, into the entire space for a suitable time interval t > 0. 

5. Some Examples 

Now let us consider several elementary solutions of the characteristic equation fl2.13p : more com- 
plicated cases may include special functions, like Bessel, hypergeometric or elliptic functions pL], 
[20] . [2T] . and [27J. Among important elementary cases of our general expressions for the Green 
function (IZTHll - ( 12^221) are the following: 

For the traditional diffusion equation 

du u 

— - = a—-T, a = constant > (5.1) 

ot ox^ 

the heat kernel is 

Equation (14. ip gives the steady solution uq =constant for all times t > 0. See [3] and references 
therein for a detailed investigation of the classical one-dimensional heat equation. 

The diffusion-type equation 

du d'^u „ , ^, 

^ = «^ + /^«' (5-3) 

where a > and / are constants (see [7], [8], [9], [10], [11], p] and references therein regarding to 
similar cases of the Schrodinger equation), has the the characteristic function of the form /i = 2at. 
The heat kernel is 

K (x, y, t) = -y^= exp | _i:LZlL j exp ( - (x + y) t + —t^ ] (5.4) 
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provided t > 0. Evolution of the uniform initial data u {x, 0) = Uq =constant is given by 

M(x,t)=Moe^"*+'^^'*'/^ (5.5) 

The initial value problem for the following diffusion-type equation with variable coefficients 

dzL lOu \ ( dzi\ 

— = a ( — ^ -x^u\ +LU i cosh ((2a - 1) t) xu + sinh ((2a - 1) t) ^ 1 , (5.6) 

where a > and uj are two constants, was solved in [15] by using the eigenfunction expansion 
method and a connection with the representations of the Heisenberg-Weyl group N (3) . Here we 
apply a different approach. The solution of the characteristic equation 

fx" - 4aV = (5.7) 
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is /i = sinh (2at) and the corresponding heat kernel is given by 

1 / {x'' + y^)cosh{2at)-2xy \ 

* '"■ "' *' = ^2.AnH2at) "''•' [ ^^^i±W) ) <'■'' 

/ X sinh (t/2)+w sinh ((2a- 1/2) t) , ft 

X exp 2c<; ^ smh ' 

V sinh (2at) 

2Cosh(2at)^^^^4/t 



Indeed, by ( l2Tfll -( l2A9l) 



X exp —2uj sinh - 

^' sinh(2at) V2 

X exp ( ^ ( t - 2 sinht + - sinh (2t) ) ) , t > 0. 



a = ^ = _^£!M2at^ g ^ 1 .5 9^ 

' 2sinh(2at)' sinh (2at) ' ^ ' ^ 



In this case 



2a 



= uj (cosh ((2a - 1) t) sinh (2at) - sinh ((2a - 1) t) cosh (2at)) 
= uo sinh t 
and equation (I2.20p gives 



By (I22ID 



cosht-1 sinh^ (t/2) 

d) = UJ — = 2uj -. (5.10) 

sinh (2at) sinh (2at) 



1-cosht 
sinh (2at) cosh (2at) 



, 1 - cosh r , f cosh ((2a - 1) r) , 
+2auj / ^- dr + uj ^V^ ^^ dr 



cosh^ (2ar) ' Jo cosh (2ar) 



where the integration by parts gives 



1 - cosh T , , , , sinh (2at) f sinh (2ar) 

dT = {I- cosht) -4 { + / -4 7 sn 

cosh (2ar) cosh (2at) Jo cosh (2ar) 



Thus 

cosh (2at) /"* sinh (2ar) sinh r + cosh ((2a — 1) t 



e = io {1 — cosh t) — -— +00 r dr 

sinh (2at) Jo cosh (2ar) 

and an elementary identity 

sinh (2at) sinht + cosh ((2a — 1) t) = cosh (2at) cosht (5-12) 

leads to an integral evaluation. Two other identities 

cosh (2at) cosht — sinh (2at) sinht = cosh ((2a — 1) t) , (5.13) 

cosh (2at) - cosh ((2a - 1) t) = 2 sinh (t/2) sinh ((2a - 1/2) t) (5.14) 

rGsiilij m 

cosh (2at) -cosh ((2a- l)t) sinh (t/2) sinh ((2a - 1/2) t) 

e = UJ = 2uj ^ . (5.15) 

sinh(2at) smh (2at) ^ ^ 
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In a similar fashion, 



9 , 4 / , N cosh (2at) 1 9 / , 1 



K = -2uj^ sinh* (t/2) — -^ f + -cjM t - 2 sinht + - sinh (2t) , (5.16) 

^ ' ' sinh (2at) 2 V 2^7' ^ ^ 

and equation (I5.8P is derived. In the limit tu ^ this kernel gives also a familiar expression in 
statistical mechanics for the density matrix for a system consisting of a simple harmonic oscillator 

m. 



The case a = 1/2 corresponds to the equation 



du 1 fd'^U 2 1, /c; 1^7^ 

a; M + a; xw l^-l' J 



at 2 V5x2 
and the heat kernel (15.81) is simplified to the form 

e-^^/^ / {{x-ujf + {y-ujf)cosht-2{x-uj){y-uj) \ 

K(x,y,t) = = exp — ^ ; , 5.18 

^ ^ ^ V27rsinht ^l 2sinht '' ^ ^ 



when t > 0. A similar diffusion- type equation 

+ X^M \ + U! xu (5.19) 



du 1 f d^u 



dt 2 \dx^ 
can be solved with the aid of the kernel 

V27rsint ^ 2smt J 

provided < t < 7i/2. We leave the details to the reader. 

Following to the case of exactly solvable time-dependent Schrodinger equation found in p!7], we 
consider the diffusion-type equation of the form 

-— = cosh^ t —— + sinh^ t x'^u + - sinh 2t 2x^— + u] . (5.21) 

ot ox^ 2 \ ox J 

The corresponding characteristic equation 

/x" - 2 tanh t /i' + 2/i = (5.22) 

has two linearly independent solutions 

^i = cos t sinh t + sin t cosh t, (5.23) 

fi2 = sin t sinh t — cos t cosh t (5-24) 

with the Wronskian W (ni, 1x2) = 2cosh^t, and the first one satisfies the initial conditions (I2.15p . 
The heat kernel is 

K{x,y,t) = = (5.25) 

^/2n (cos t sinh t + sin t cosh t) 

iy"^ — x"^) sin t sinh t + 2xy — (x^ + y^) cos t cosh t 
X exp ' 



2 (cos t sinh t + sin t cosh t) 

provided < t < Ti ^ 0.9375520344, where Ti is the first positive root of the transcendental 
equation tanht = cott. Then 7 (t) < and the integral (I2.26P converges for suitable initial data. 
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A similar diffusion-type equation 

du n d'^u o n 1 f du \ , ^ 

—- = cos^ t — - + sm^ t x^u sm 2t 2x—- + u (5.26) 

at ox^ 2 \ ox J 

has the characteristic equation of the form 

/i" + 2tant /i'-2/i = (5.27) 

with the same solution (15.231) . It appeared in [IT] and [6] for a special case of the Schrodinger 
equation. The corresponding heat kernel has the same form (I5.25P but with x and y interchanged: 

K{x,y,t) = = (5.28) 

^/2-K (cos t sinh t + sin t cosh t) 

(x"^ — y"^) sin t sinh t + 2xy — ix'^ + y'^) cos t cosh t 
X exp ' 



2 (cos t sinh t + sin t cosh t) 

provided < t < T2 ~ 2.347045566, where T2 is the first positive root of the transcendental equation 
tanht = — cott. We leave the details for the reader. 

6. Solution of the Non-Homogeneous Equation 
A diffusion-type equation of the form 

^-g(t)^« = F, (6.1) 

where Q stands for the second order linear differential operator in the right hand side of equation 
(12. ip and F = F [t, x, u) , can be rewritten formally as an integral equation (the Duhamel principle; 
see [1], [S], [13], [II], [21], [2S] and references therein) 

ti (x, t) = if (t, 0) M (x, 0) + / H{t,s)F{s,x,u) ds. (6.2) 

Jo 

Operator H{t,s) is given by (I2.28p . When F does not depend on u, one gets a solution of the 
nonhomogeneous equation (16. ip . 

Indeed, a formal differentiation gives 

du d d f^ 

-g^ = g-^H (t, 0)u{x,0) + -j H (t, s) F {s, X, u) ds, (6.3) 



where 

d_ <" 

dt 
and we assume that H (t, t) is the identity operator. Also 



/ H{t,s) F{s,x,u) ds = H{t,t) F{t,x,u)+ / —H{t,s)F{s,x,u) ds (6.4) 

Jo Jo (Ji 



and 



Q(t)u = Q(t)if (t,0)M(x,0)+ / Q{t)H{t,s)F{s,x,u) ds (6.5) 

Jo 



^-Q(i))^ = (^^-Q{t)]H{t,0)u{x,0) + F (6.6) 
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+ / (^-g(t)W(t,s)F(s,x,n) ds, 



where 

d 



^^ Q{t)\H{t,s)=0, 0<s<t (6.7) 

by construction of the operator H {t, s) in (12.281) . This completes our formal proof. A rigorous proof 
will be given elsewhere. 
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